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Abstract 

In this chapter the recent theoretical work on phase transition in imbalanced fermion superfluids 
is reviewed. The imbalanced systems are those in which the two fermionic species candidate to 
form pairing have different Fermi surfaces or densities. We consider systems subjected to weak 
interactions. In this scenario two distinct phase transitions are predicted to occur. A thermody- 
namical phase transition, induced by the temperature (T), and a quantum phase transition as a 
function of the increasing chemical potentials asymmetry, that takes place at zero temperature. 
We also briefly discuss some recent experimental work at non-zero T with imbalanced Fermi gases 
in cold atomic traps. 
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I. INTRODUCTION 



It is well known that fermions with opposite spins near a common Fermi surface suffer the 
Cooper pairing instability at arbitrarily weak coupling, below a certain critical temperature, 
leading to the phenomenon of superfluidity. When the two fermion species have the same 
density, the ground state is described by the successfully Bardeen-Cooper-Shrieffer (BCS) [l[ 
theory of superconductivity A novelty in the pairing mechanism is brought about by an 
imbalance in the number densities of the two species 1 : certainly there will be fermions 
without a partner. 

For this asymmetrical scenario various exotic phases have been suggested, such as the first 
proposal 2 by Sarma [H, the Larkin, Ovchinnikov, Fulde and Ferrel (LOFF)-phase j| (where 
the gap parameter breaks translational invariance), deformed Fermi surfaces jlfj, [ll| (in 
which the Fermi-surfaces of the two species are deformed into ellipsoidal form at zero total 
momentum of Cooper pairs), the breached pair superfluid phase (BP) [l2|, E3 (composed 
by a homogeneous mixture of normal and superfluid phases), and phase separation in real 
space [T5I . Eq l (defined as an inhomogeneous mixed phase formed by normal and superfluid 
components). Alternatives beyond mean field have also been presented, such as the (induced) 
P-wave superfluidity fl7^ . 

The very recent experimental demonstration of superfluidity through the observation of 
vortices and phase separation with controlled population imbalance in atomic gases 0, 
2ol ] have promoted a great interest in the field, including prospects from the theoretical point 
of view 21, 22I 23|, 2J, 2^, 26|, 27, 28|. In these experiments, the strength of the interactions 
can be tuned by the use of the Feshbach resonance, varying the external magnetic field. In 
this way, the crossover from BCS superfluidity to the Bose-Einstein condensation (BEC) 
can be accessed. For a more complete description on recent experimental work on various 
novel superfluid phases in fermion systems, see the chapter Realization, characterization, 
and detection of novel superfluid phases with pairing between unbalanced fermion species in 
this volume, Ref. [29l |. 

In passing we note that the counterpart for ordering in the quantum chromodynamics 
(QCD) context (for the case where the quarks have approximately the same Fermi surface) at 
high baryon density, is the color superconductivity jflU 0, H3| • For high density asymmetric 
quark matter the analog of LOFF 3 state leads to crystalline color superconductivity jjjjj. In 
QCD with two light flavors, the ground state is the two-flavor color superconductor (2SC) 
[34| . The pairing between u and d quarks under the condition of charge neutrality and (5- 
equilibrium at intermediate baryon densities, denoted gapless color superconductivity, has 
been studied in Refs. jlH |36| . The quark (and nuclear) pairing will be discussed in other 
chapters of this issue. 



1 These different species could be two fermionic atoms ( 40 K or 6 Li) or hyperfine states of the same atom 

The nucleus of neutron stars could also have the basic ingredients for asymmetric pairing: quarks 
with different flavors |EIE0- 

2 Sarma considered a model exhibiting asymmetry between the (fixed) chemical potentials of two particle 
species at zero and finite temperature |8j. 

3 A review of the theoretical approach and phenomenological applications of the LOFF state in Condensed 
Matter and QCD, can be found in Ref. [31|. An analysis of the competition between the BP and LOFF 
pairing mechanisms in asymmetric fermion superfluids is shown in Ref. |.'S2| . 
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II. THE MODEL 



To begin with, let us consider a nonrelativistic dilute (i.e., the particles interact through a 
short-range attractive interaction) cold fermionic system, described by the following Hamil- 
tonian 

'K = H ^n a = ^ e k a l a k + e b k b{b k - g ^ a \>b ] - k ,b„ k a kl (1) 

k,ct k k,k' 

where a\, a k are the creation and annihilation operators for the a particles (and the same for 
the b particles) and e k are their dispersion relation, defined by e k = £, k — /i a , with = ^— 
and fi a being the chemical potential of the (non-interacting) a-specie, a = a, b. To reflect 
an attractive (s-wave) interaction between particles a and b we take — g < 0. 



1. The Thermodynamic Potential 

We now derive the thermodynamic potential at fixed chemical potentials 4 , and the finite 
temperature gap equation for an asymmetrical fermion system, in order to determine the 
critical temperature. We follow the usual derivation of the textbooks [371 ]. however extending 
the analysis for the imbalanced systems we are investigating. The thermodynamic potential 
is given by 

n = E-TS, (2) 

where E is the internal energy and S the entropy. Let us define f k as the probability of an 
a particle with momentum k is excited, and similarly g k as the probability of a b particle 
with momentum — k is excited. One can write the following possible probabilities for the 
particles states: The probability that a given pair of k is unexcited is 

P fc (0) = (l-/ fc )(l -(/*)• (3) 
The probability that one of the states (k, for instance) is excited, and the other is not, is 

P k (2) = f k (l-g k ). (4) 
Now the probability that the state — k is excited, and the other is not, is 

^(3) = (1 - fk)9k- (5) 
And finally, the probability that both states are excited is 

Pk(l) = fk9k- (6) 

The entropy is defined as 

S = — ^^probability of state i x In [probability of state i] = — Pi ln(Pj), (7) 



4 The situation where the chemical potentials are kept fixed can find place, for instance, in a gas of fermionic 
atoms connected to reservoirs of species a and b, so the number densities are allowed to change in the 
system. 
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where we have set the Boltzmann constant equal to one. It is left as an exercise to the reader 
to show that S for an asymmetrical fermion gas is found to be 



S = ~ E {/* ln ^) + - A) ln ( X - A) + ^ ln (^) + ( X " -9k)}- (8) 

The BCS ground state, which describes a superposition of empty and occupied (paired) 
states, is given by [l[ 

\BCS) = H[u k + v k albl k ] |0>, (9) 

k 

where the arbitrary complex (a priori) coefficients u k and v k are to be determined by a 
variational calculation. They are subjected to the normalization, \uk\ 2 + \vk\ 2 = 1, and 
spin singlet, u k = U- k , v k = V- k , conditions. At zero temperature, the internal or ground 
state energy E is simply the expectation value of the Hamiltonian, E =< BCS|CK|BCS >. 
At finite temperature the energy has to take into account the excitations probabilities, 

E = Efc E i p i & Then we find 

E = E {3K 1 " /* " M + /*] + 4K1 - /* - 9kH + g k ]} (10) 

k 

-9^2 u k'Vk'U k v k (l - fk- gk)(l - fk> - 9k>)- 



k,k' 



Plugging Eqs. (JHJ) and (|10|) in Eq. (J2J), and taking the minimizations 



we find, respectively, 



0, (11) 
oh 

i°=o, 

8 9k 

|°=0, 



A = l/(e" e 8 + 1), (12) 
^ = l/(e^ + l), (13) 



2 _ i A , 4 
k 



where 



2V 1 + tJ' (M) 
E? = ±ej; + E k , (15) 
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are the quasiparticle excitations, with E k = J e^ 2 + A 2 (T) and = efc ^ €fc . In terms of f k , 
g k and u k the thermodynamic potential becomes 

n=— + J2[4-E k - Tln(e-^ + 1) - Tln(e"^ + 1)] . (16) 
9 k 

In the definition of 8.^ b we have also defined 

A{T)=g^2u k v k (l-f k -g k ). (17) 

k 

Since v k = 1 — u\, then = and the gap equation can be written as 

l =9Y,w {1 ~ h ~ 9k) - (18) 

k 

We note that the equation above can also be obtained by |£ = 0. 



III. PHASE TRANSITIONS 

A. Phase transitions in fully gapped systems 

The critical temperature T c is, by definition, the temperature at which A = 0. Then 
Eq. (fTHj) becomes 

(19) 

After the momentum integration, Eq. (|19|) can be written as [3£ 

_L__ ln ( A ^)._I y(a) , (20 ) 

where p(0) = is the density of states at the Fermi level, with kp = y/2Mfi being the 
"average" Fermi surface having fi = fi a + /x b , and M = ™ a ™^ b is the reduced mass. We 

also introduced o = . = as a dimensionless parameter reflecting the mass or chemical 

potential asymmetry, and ^(a) = ^(| + —) + ^(| — —), with ^ being the digamma function, 

defined as ^(z) = where z is a complex number with a positive real component, T is the 

gamma function, and V is the derivative of the gamma function and a= — ^b^b-maHa 

° ' ° 2 ' 2 in a +irif, 

Eq. (|2()jl can be put in the form 

T c =^e-i*M (21) 

where a c = yV an d = 2c<jd e^ 1 ^^ 9 is the BCS gap parameter in the weak coupling limit, 
p(0)g << 1. The critical temperature we have obtained refers to the situation where all the 
fermions are gapped 5 ; they have the same Fermi surface in spite of the masses and chemical 



By construction the a and b particles in the BCS state have the same number densities, see Eq. @- 



5 



0.60 



0.55 



0.50 



0.45 



0.40 



0.35 



0.30 




FIG. 1: T c /Aq as a function of rrib/m a for a system constrained to m a (i a = rribUb- 



potentials asymmetry. The consequence is that P F = P F (where P F ' = ^2m a ^^ a ,b) implies 
a c = 0, resulting in 



T c {P a F 



2a— A , 



7T 



(22) 



where we have used that 3^(0) = —27—4 ln(2), where 7 is the Euler's constant. An important 
feature of Eq. (J2~2"j) is that, although Eq. (fT§|) requires regularization, the regulator depen- 
dence cancels from the result (}2"Tj) [381 ] . In the symmetric limit, namely m a = nib (or /i a = fit,) 
the standard BCS relation T c /A = — [3, H3] is recovered. We can observe from Eq. (|2*2*|) 
that the critical temperature for the system constrained to P F = P F (or m a /i a = m b fi b ) goes 
with 2y/^ ^Ao for greater than m a and approaches zero for mj >> m a . This show 

that the pair formation is disfavored for very large mass asymmetry, even in systems were 
the Fermi surfaces match. In Fig. we show the ratio T c /A as a function of mb/m a . 
As one can see, T c /Ao is a smooth function of the mass asymmetry, and goes to zero for 
m b /m a -> 00. 



B. Phase transitions in gapless systems 



Depending on the relative magnitudes of the particles Fermi surfaces and masses, the 



quasiparticle excitations 8.^ b can be negative. If we choose P F > P F and m& > m„. onl y E j 
will cross zero. From the equation ele b k = —A 2 we determine the roots of E b k 



16m a m fe A 2 ] 1/2 , 



(23) 



where 5Pp = Pp* ± Pp 2 - The negativity of <L b k between k\ and k 2 means that the corre- 
sponding states (the b particles) are singly occupied. In Fig. (J2J) we show the quasiparticle 
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excitations (QPE) as a function of the momentum k. For some values of A, £^ may be 
negative for momenta k\ < k < /c 2 . 

We can observe from Eq. (J23j) that A has a critical value 

\fiP~ 2 \ 

A c = /111 . (24) 

For A > A c , are not real and £^' 6 never crosses zero. This corresponds to the standard 
BCS with pairing for all k. The situation where A < A c , named "Sarma phase", was first 
pointed out in Ref. 0. 

The thermodynamic potential in the Sarma phase is obtained when we find a state 
which minimizes the internal energy. The smallest energy is reached when the modes with 
negative 8, k ,b are filled and the remaining modes are left empty. More precisely, the ground 
state |$) satisfies 0,Gi| 

a k ,b k \V) = if ££>0, 



,t i.t 



b\m = if £ b fc <0. (25) 



This state can be written in terms of the a\ and b k operators and the vacuum state |0) as 



i*) = n [' 



k-2 



u k + v k a\bl k \ J6t| ). (26) 



k<k± fei 
k>k 2 



The state above corresponds to having BCS pairing in the modes k where £^ > and a 
state filled with particles b in the modes where E\ < 0. Using this state in the computation 
of the entropy and the internal energy, the thermodynamic potential of the Sarma phase 
turns out to be 



Q(g,T,A< A c ) = Q S (T) (27) 
= —+J2[4-E k -T\n(e-^ + 1) - Tln(e"^ + 1] 

^ k<k! 
k>k 2 

k 2 



+ J2[e b k -THe^ + l) 



Since there are gapless states in Eq. (J2"7j) one can not define a critical temperature as 
performed in the previous situation, i.e., in the fully gapped system. To find the critical 
temperature (T*) in this case, it is necessary a comparison between Q S (T) and the normal 
free energy, Q^fA = 0,T), at a given and fixed asymmetry 5ji < 5fi c 6 , for increasing 
temperature j41(. We then define T* for gapless systems as the temperature at which 
n(SfM, T = T*, A = A ) = n(6fi, T = T*, A = 0). Solving this equality for T* we obtain the 
transition temperature. 



6 The prediction for the break down of the fermionic superfluidity is Sfi c = Mb 2 A '° = Ao |l5l Il6l l57j . This 
picture has been confirmed qualitatively experimentally . 
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FIG. 2: Dispersion relation for the quasi-particles a and (3 showing a region where £^ is negative 
for m b = 7m a , P b F = 1.45 PJ, and A = 0.29A , obtained by Eq. fS%). Solid curve corresponds to 
£^ and dashed curve corresponds to 



We point out here that the effect of the temperature in the thermodynamic potential of 
Eq. (|2"7|) is to induce a transmutation from a stable phase at Q(Sfi,T < T*,A = A ) to an 
also stable phase Q{Sfi,T > T*,A = 0) HJ 0- The nature of this transition, however, 



needs further investigation and will be presented elsewhere |41| . 

C. The quantum phase transition 

The zero temperature limit of Eq. (|27|) yields 



A s 



k 2 



Q(g, T = 0, A < Ac) = Q S (T = 0) = — + £ [e+ - E k ] + ^ e 



(28) 



k<ki 
k>k 2 



With 8n if °' ) = (remembering that the partial derivative also hits fci and &2 in the limits 



0A 



of the integrals) we obtain the gap equation 



d 3 k 



2 ><fci (2tt) 3 



4 + A 2 



We can find the gap in the Sarma phase through the identity [1, HI 3, 



M 
1ix\a\ 



d 3 k 



d 3 k 



+ A^ >/ fc>fc 2 v y \/e7 + A 2 



(29) 



(30) 



For small values of the gaps (A , A << /z a , the integrals can be approximated and it is 
found that 

A 2 eKh) 



A 2 , et(k 2 ) 



(31) 
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FIG. 3: Thermodynamic potential for different values of P F and P F (constant kp). The top curve 
corresponds to Pp = Pp and the lower curves correspond to increasing values of 



\pb 

' F 



pa I 



which has the solution 



\ 



Ao(Jpl-Ao). (32) 
1 2 v /m a m b 



It is ease to verify that the Fermi surfaces asymmetry and the gap in the Sarma phase are 
restricted to: 

2y/m a m b A < \SP F 2 \ < Ay/m a m b A (33) 
< A s < A . 

In Fig. (J3J) we show the thermodynamic potential as a function of A for different values 
of P F and Pp, keeping the combination kp/M = P F 2 /m a + Pp 2 /m b fixed, computed from 
a numerical evaluation of Eq. (|28|) . As we can see, there exist a special combination of Pp 
and P F for which Q S (T = 0) has double minima. Then, given a Pp, we want to know the 
correspondent value of P F that satisfies this requirement. We note from Eq. (J23|) that when 

A = 0, fci = P F and k 2 = P b F , whereas for A = A , h = k 2 = \J 1 p f 2 + p f . . Thus, the 
condition to find the relation between P F and P F is 

n s (T = 0, A = 0) = n s (T = 0, A = A ) ->■ (34) 

Note that the l.h.s. is just the free-energy of the unpaired a and b particles and the r.h.s. is 
the standard BCS free energy. Solving both sides we get 

where k F = [M(Pf/m a + Pp 2 /m b )} 1/2 , and A = ^e^^ is the standard BCS gap pa- 
rameter [3]. The solution of the equation above we define as P FD mi where DM stands 
for double minima. Given these fixed values of P F and P F (or \x b and ji a ) found above, the 
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system could either be found in the normal, BCS or in a mixed phase, since those states 
would have the same minimal energy. When the free-energy F(P FDM ) has its minimum at 
A , after an increase in P F the gap jumps from A to 0, characterizing a first order quantum 
phase transition from the superfluid to the normal phase, as has been found theoretically in 
several physical situations. See, for instance, Refs. [42|, |43|, |46| . Summarizing, we find that 
at fixed chemical potentials the Sarma gap corresponds to an extremum of the free-energy, 
it never represents its minimum. 

For any curve having momenta P F > Pp (and fixed in each curve of Fig. Q) until PpDMi 

both particle species have the same density n a = n b = — m d ^~^ = ^ + R^f 
For P F > P F 

DM-* there is absorption of particles from the reservoirs so the asymmetry in 
the number densities gets its maximum value, and the system goes to the normal state with 

number densities n a = and n b = -^y. 

The mixed phase (MP) or heterogeneous composition is formed by a normal and a super- 
fluid components, and (differently from the two ho mogene ous ground states found above) 
accommodates n a and n& particle densities in a trap 
stable, with no surface energy cost at weak coupling 
in the MP is under investigation jiT 



THE. The MP has been found to be 
44j . The issue of the phase transitions 



IV. BRIEF REVIEW ON SOME ASPECTS OF RECENT EXPERIMENTAL 
WORK ON FERMI GASES AT FINITE T IN ATOMIC TRAPS 



For unitary Fermi gases, current experiments [421, |4ja, |4jJ, |£0j produce temperatures down 
to about 0.05Tp, where Tp is what the Fermi temperature would be for a noninteracting 
gas with the same number of atoms and in the same trap conditions. Typically Tp is 
of order fiK. However, a weakly interacting Fermi gas requires much lower T to achieve 
superfluidity. For the conditions of these experiments, the mean field approximation with 
an interaction energy proportional to the scattering length is not valid. However, the mean 
field approximation with a unitary limit appears approximately valid, furnishing a good 
agreement with predictions of the collective frequencies, and a very good agreement on the 
transition temperature |5l| . 

Recently, measurements of the T-dependent momentum distribution of a trapped Fermi 
gas consisting of an equal mixture of the two lowest spin states of 40 K in the BCS-BEC 
crossover regime have been presented jH^ |. The results show the existence of a competition 
between the T dependence of the fermionic excitation gap and thermal broadening, leading to 
non-monotonic behavior in the T dependence of the momentum profiles. Semi-quantitative 
agreement between theory and experiment using a simple mean-field theory has been found. 

Thus, even when the measurements are done in strongly interacting Fermi gases, mean 
field theory had qualitatively explained the behavior of these systems, and we expect that 
our weak coupling mean field BCS results should also be valid, at least, qualitatively. 

More recently, direct normal-to-superfluid phase transition has been observed in a 
strongly interacting Fermi gas with unequal mixtures of the two spin components jH^. 
Both the thermodynamical and quantum phase transitions have been detected. A quantum 
phase transition at 5 ~ 70% has been observed, where 5 = is the imbalanced param- 



eter, with iVj(iVj) as the majority (minority) atom number. The authors of Ref. |53[ found 
that the critical temperature will in general depend on the population imbalance. This is a 
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(expected) guide for theoretical investigations |4 



V. SUMMARY 

In this chapter we have discussed phase transitions in cold fermionic gases composed 
by two particle species whose Fermi surfaces or densities do not match. We have seen 
that depending on the relative difference between the Fermi surfaces, two distinct phase 
transitions are predicted to occur. 

Fully gapped asymmetrical systems undergo a second order phase transition, driven by 
the temperature, obeying ^P- = J ^^- « 1.76 77-, where a = ^ m °^ rib = appears to be 

an universal constant |55| . Since experiments are being set up to study pairing between 
fermions of unequal mass (for example 40 K and 6 Li), we expect that this expression can be 
verified soon |56j]. 

A system with different Fermi surfaces is found to be in the BCS state while the asymme- 
try between the chemical potentials is smaller than the critical difference Sfi c . A quantum 
first order phase transition from superfluid to normal phase, at zero temperature, happens 
as a function of the increasing chemical potentials asymmetry. 

As systematic studies at non-zero temperature in imbalanced ultracold systems are just 
starting, we hope that the work presented in this chapter could stimulate new investigations 
in this interesting field. 
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